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ABSTRACT 

We present a theoretical model for the thermal X-ray emission properties and cool- 
ing behaviors of isolated pulsars. It is assumed that pulsars are solid quark stars. 
Using Debye elastic medium theory, we calculate the heat capacity of such a quark 
star and show that the residual thermal energy cannot sustain the observed thermal 
X-ray luminosities seen in typical isolated X-ray pulsars. We conclude that another 
heating mechanism must be in operation if the pulsar is in fact a solid quark star. 
Two possible alternative heating mechanisms are explored. First, for pulsars with lit- 
tle magnetospheric activity, accretion from the interstellar medium or from material 
in the associated supernova remnant may power the observed thermal emission. A 
disk-accretion rate M ~1 % of the Eddington rate, with an accretion onto the stellar 
surface at a rate of ~ 0.1 %M, could explain the observed emission properties of dim 
isolated neutron stars as well as central compact objects. Second, for pulsars with 
significant magnetospheric activity, the spindown of the pulsar itself may be providing 
the source of thermal energy via the reversing plasma current flow. A phcnomenolog- 
ical study between pulsars' bolometric X-ray luminosities and their spin energy loss 
rates presents the probable existence of a 1/2-law or a linear law, i.e. L^ ol oc E 1 / 2 or 
L^ ol oc E. This result together with the thermal properties of solid quark stars allows 
us to calculate the thermal evolution of such stars. Cooling curves are calculated and 
compared with the thermal emission properties measured from 17 active X-ray pul- 
sars. It is shown that the observed properties of these sources are consistent with the 
solid quark star pulsar model. 
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1 INTRODUCTION 

The study of cold quark matter phases has been a hot 
topic of research in recent years. In an astrophysical con- 
text, quark matter may be used to explain the measureable 
properties of pulsar-like compact stars. 

It has recently been proposed that real i stic quark mat- 
ter could exist in a solid state JXul 120031: iHorvathl I200I 
IOwenll2005l : iMannarelli, Raiagopal fc Sharmal 120071 '), either 
as a super-solid or as a normal solid l|Xul 12003 ). The basic 
conjecture of normal solid quark matter is that de-confined 
quarks tend to form quark-clusters when the temperature 
and density are relatively low. Below some critical tempera- 
ture, these clusters could be in periodic lattices immersed in 
a degenerate, extremely relativistic, electron gas. Note that 
even though quark matter is usually described as weakly 
coupled, the interaction between quark s and gluons i n a 
quark-gluon plasma is still very strong (|Shurvakl 120091 ). It 
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is this strong coupling that could cause the quarks to clus- 
ter and form a solid-like material. 

In this paper, we argue that comtemporary observa- 
tions of thermal X-ray emitting pulsars are consistent with 
the assumption that these sources are in fact Solid Quark 
Stars (SQSs). As it turns out, the heat capacity of solid 
quark matter can be easily estimated using Debye elastic 
medium theory. This is done in section §2. It is also shown 
that the internal energy of a SQS would be radiated almost 
instantaneously after its birth. Hence, if pulsars are SQSs, 
the observed thermal X-ray radiation must be powered by 
some other sources of energy. Two possible energy sources 
are proposed. The observing consequences of these models 
are discussed in 83. Our conclusion are discussed in 84. 



2 THE MODEL 

2.1 Cooling stages of a quark star 

The cooling process of a quark star can be quite complicated 
to model when starting from the birth of the star. Theo- 
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Figure 1. Possible cooling stages of a quark star. Stage 1: individual quark phase. The temperature is high (> 10 MeV) when the star 
is born, and the state of the star could be the fluid of individual quarks. Stage 2: individual quark cluster phase. As the temperature 
decreases, individual quarks tend to form quark clusters because of the strong coupling between them. The state of the star could then 
be the fluid of quark clusters. Stage 3: solid quark star phase. As the temperature drops to the melting temperature, the fluid of quark 
clusters tend to solidify to form periodic lattice structure, such as bcc structure. 



retical uncertainties make it difficult to predict the exact 
temperatures where phase transitions occur. For illustrative 
purposes, we present the following scenario where a SQS's 
cooling takes place in approximately three stages (See Fig. 
1). The first stage occurs just after the birth of the quark 
star if its initial temperature is much higher than ~ 10 11 K 
(10 MeV). The emission of neutrinos and photons will lead 
to fast cooling. Hence, the star quickly enters the second 
stage where de-confined quarks begin to form quark clus- 
ters. As the temperature drops further, the fluid solidifies 
and the star is said to enter the third cooling stage. Here, 
quark clusters form periodic lattice structures (e.g. the bcc 
structure). 

It is uncertain exactly how long the star will spend in 
each phase, or even if it will ever enter all three states. If 
the initial temperature is only just around 10 11 K or lower, 
Stage 1 may be short lived or even non-existent. The quark- 
gluon plasma could be strongly coupled at birth and quark- 
clusters could be present immediately after formation. It 
is even possible that the melting temperature of solid quark 
matter could be 10 11 K or higher. Hence, a SQS would enter 
Stage 3 immediately. 

The mechanisms for the emission of thermal photons 
for each specific stage could be quite different. Hot quark 
stars (in Stages 1 or 2) could be a good radiator for ther- 
mal equilibrium photons with energies more than ~ 20 MeV 
(|Alcock. Farhi fc Olintol I1986T ). Meanwhile, the intense re- 
lease of thermal energy could stimulate the generation and 
radiation of electron-positron pairs from the hot bare quark 
surfaces; nonthermal photons generated by electron-electron 
bremsstrahlung processes could also be an emitting compo- 
nent though the amount could be small (|Usovll200lh . When 
the stars cool down and become solidified (in Stage 3), ther- 
mal equilibrium photons (usually soft X-ray) could also be 
emitted from bare SQSs' surfaces as free electrons tran- 
sit in the levels of the energy bands of solid quark matter 
|Xull2003h . The electron-phonon interaction as well as the 



interaction between electrons themselves could result in a 
metal-like spectrum. IZhang. Xu fc Zhand l|2004l ) have ana- 
lyzed this spectrum phenomenologically, but they did not 
conclude that there are significant differences between the 
metal-like spectrum and black body spectrum. Hence, black 
body thermal spectrum could be a good approximation for 
SQSs. 

This work concentrates its attention on the late third 
stage with surface temperature varies between ~ 10 7 and 
~ 10 K, which is the range of thermal components of iso- 
lated pulsars that telescopes have detected. The neutrino 
cooling could be negligible for a SQS in the late third stage, 
since the temperature then is low enough and the pho - 
ton luminosity may dominate |Yakovlev fc Haensell 120030 . 
This situation is quite different from a normal neutron star, 
which i s surrounded by a nearly adiabatic blank et (i.e., the 
crust; iGudmundsson. Pethick fc Epstein! 1 19831 ). Hence, in 
this stage, the stars cool down by losing thermal photons, 
and the sources of the energy that enable the thermal emis- 
sion are firstly the residual thermal energy of the stars and 
secondly the energy input caused by the processes of stellar 
heating. The energy relation could then be described as 



dU 
dt 



+ Lsh = Lb 



(1) 



where Lboi is the stellar bolometric X-ray luminosity, 
—dll/dt is the release rate of the stellar residual thermal 
energy, and Lsh is the luminosity of stellar heating. 



2.2 Residual thermal energy of a SQS 

Provided that the volume of a SQS is a constant, then the 
stellar residual thermal energy Usqs could only be the func- 
tion of the star's temperature T s , 



Usqs{T s , 



C v dT Sl 



(2) 
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where T s is the value in the star's local reference frame, and 
C v is the heat capacity of the star. Following Debye elas- 
tic medium theory, the characteristic of the heat capacity of 
solid quark matter could be evaluated by Debye tempera- 
ture, 



6 D = hu) D /kB, 



(3) 



where h is reduced Planck constant, and fca is Boltzmann 
constant, uoo is Debye cut-off frequency (i.e. the maximum 
frequency of the wave that could propagate in a medium), 
which equals Debye wave number kn = (67r 2 n c ) 1//3 (n c is the 
number density of classical particles, or quark clusters for 
solid quark matter especially) multiply the average sound 
speed in the medium, i.e. u>d ~ koc 3 . For a SQS, the average 
sound speed could be the light speed approximately. A linear 
equation of state, extended to be used for a quark star, notes 
that the pressure p ~ pc 2 , where p is the mass density of a 
quark star, and c is the speed of light, so an estimate could be 

V d P/ d P 



c a = y/ dp/dp ~ \JpIp ~ V 'pc 2 ' I ' p — c. n c = 3eno/A, where 
e denotes the baryon number density of solid quark matter 
in the unit of no • no is the baryon number density of normal 
nuclear matter and equals 0.17 fm -3 . We consider e = 3 — 5 
could be the representative values for a SQS. Noting that in 
the following calculation in this paper, we will adopt e = 3 
as a typical value, since the variation of e between 3-5 could 
not lead to the magnitudes of the results varying in orders. 
A is the number of valence quarks in a quark cluster. We 
may expect that A could be in the order of 10, since A = 1 8 
if quark-a-like clusters are formed l|Michellll988l ; ,Xu,l2003h , 
and A could even be conjectured to be in the order of ^ 10 2 . 
Debye temperature, 9d, of a SQS could then be of order 
~ 10 12 K, even higher than the temperature when a SQS 
is born. Hence, the heat capacity in the low-temperature 
limit, or the temperature-cube law, is applicable for the third 
stage, i.e. 



127T 4 



7 ( T B v3 



(4) 



where c v is the heat capacity per classical parti- 
cle (or quark cluster when refer s to solid quark 
matter) l |Vonsovskv fc Katsnelsonl [l989), and C v = Nc v , 
where N is the total number of clusters in a star. 

By referring to equation ([2]) , the residual thermal energy 
for a 1.4M Q -SQS could thus be U S qs ~ 10 4 T S 4 ergs, while 
the corresponding v alue for such a typical neu tron star is 
[/ NS ~ 10 30 T 2 ergs (jYakovlev fc Pethickl I200I ) . Therefore, 
the residual thermal energy of a SQS will always be much 
less than that of a normal neutron star; even when the tem- 
perature is as high as ~ 10 11 K, the residual thermal energy 
of a SQS is only one ten thousandth of that of a neutron star. 
So if the cooling process of a 1.4M©-SQS is only powered by 
its residual thermal energy, i.e. 



C v — — = 4ttR gT s + -ttR 10 T s , 9 
at 3 



(5) 



(R is the local stellar radius, a is Stefan-Boltzmann con- 
stant), then its cooling time scale from 10 11 K to an ex- 
tremely low final temperature 1 K is only ~ 0.2 s. Noting 
that as a result of the concerning of the high initial temper- 
ature here, neutrino luminosity for the SQS is also consid- 
ered as the se cond term in the right side of the equal sign of 
equation (0) (jYakovlev fc Hae nscl 2003). This could be the 



intrinsic distinction for SQSs from neutron stars. Neutron 
stars' cooling is mainly residual thermal energy powered, 
while SQSs' thermal X-ray emission and cooling processes 
could however be sustained by heating processes. 

2.3 Stellar heating 

Heating processes may play a significant role in the ther- 
mal evolution or realizing the visibility of pulsars in the soft 
X-ray band, and pulsars with different magnetospheric prop- 
erties may be undergoing different heating mechanisms. 



2.3.1 Spin origin 

Luminous nonthermal radiation and bright pulsar wind neb- 
ulae (PWNe), these observing manifestations may imply an 
active magnetosphere and strong star wind or relativistic 
particle flow could be ejecting from the poles of such a 
pulsar. A certain amount of backflow of such plasma in- 
duce the stellar heating, or an energy input could take 
place at the polar caps but disperse to the bulk of the star 
dCherig fc Rudermanfll977l ; lArons! Il98ll ; IZhang fc Hardind 
l200d ; IWang et alJll99Sh ~The heat flow, H, could thus be 



H 



Lsh - knrloT* T p -T s 



R 



(6) 



where a is Stefan-Boltzmann constant, r p w R^/QR/c and 
T p are the radius and the surface temperature of polar 
caps in the local reference frame respectively, R is the lo- 
cal stellar radius. To the simplicity, we adopt the one di- 
mensional approximation for the gradient of temperature, 
i.e. VT ~ (T p — T s )/R. Lsh, the luminosity of the return- 
current stellar heating, could be a function of the spin energy 
loss rate E, and could generally be in the form of a power 
law. However, a phenomenological study on the bolometric 
luminosity and E reveals that the power index could be 1/2 
or 1, i.e. Lsh = CE 1 ^ 2 or Lsh = r\E (see Appendix A). 

The metal-like interior of a SQS could imply that its 
thermal conductivity k could be the sum of partial compo- 
nents of phonons, electrons, and static impurities, or 



(7) 



where the subscripts 'p', 'e', and 'imp' denote the par- 
tial components mentioned above, respectively. The electron 
thermal conductivity, k b , could be the dom inant compo- 
nent for the solid metal (|Flowers &; Itohlll98lh . So phonons' 
contribution, k p , could be neglected, and, as a preliminary 
model, we omit the partial component of static impurities. 
The electron thermal conductivity could be written as 



1 1 
h 



(8) 



where the partial component contributed by the colli- 

sion between electrons, and K pe is the one contributed by the 
collision between pho nons and electrons. Fo r their analytic 
formulae, we refer to iFlowers fc Itohl ( 198ll ) with assuming 
that the results remain hold for SQSs. 



2.3.2 Accretion origin 

Some pulsars, showing steady long-term soft X-ray fluxes, 
being inert in nonthermal emission, lacking the proves on the 
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existences of the associated PWNe, (and sometimes) owning 
excesses of X-ray luminosities relative to the spin energy loss 
rates, could alternatively be inactive pulsar candidates. Due 
to their little magnetospheric manifestations, SQSs will not 
suggest the thermal X-ray radiation of these sources is of 
spin origin. 

One possible way to understand the origin of the energy 
is the accretion in the propeller regime. In this regime, a 
shell of atmosphere of matter may form surrounding the 
star. Matter closing to the inner boundary of the shell may 
interact directly with the rotating stellar magnetosphere, 
as a resu lt of which m ost of the matter will be expelled 
outward (|Lipunovlll992r ). Nevertheless, a certain fraction of 
the accreting material, described by the accretion efficiency 
?7acc, may diffuse starward and fall onto the surface of the 
star finally. In the propeller regime, the star's magnetosphere 
radius or Alfven radius r m could be on one hand larger than 
its corotation radius r co , but is on the other hand somewhat 
smaller than the light cylinder radius rt,, i.e. 



{ GMy/ 3p 2 /3 < ^ = ( 



R b Bl 



4n 



MV2GM 



f 7 < 



J'L 



cP 

2tt' 



(9) 



where G is the gravitational constant, M, R, and B p are 
the stellar mass, radius, and the magnetic strength at the 
poles, respectively. The accretion rate M could be scaled 
by Eddington accretion rate MEddj so M could denote the 
accretion rate in this unit, i.e. M = jiM^dd- During the ac- 
cretion of the material which could reach the stellar surface 
eventually, the energy of gravitation will be released. When 
the two-flavor baryonic matter impact upon the surface of 
a SQS, it will burn into the three-flavor strange quark mat- 
ter phase, and the latent heat of Ae ~ 10 MeV - ~ 100 
MeV per bary on could be released in the phase transition 
l|Madsenlll999h . The luminosity of stellar heating in this sit- 
uation could then be 

T GM ■ 77 acc ■ M 7/acc • M 

Lsh = ^ h Ae- 



R 



where m p is the mass of a proton. 



(10) 



3 OBSERVATIONS VERSUS EXPECTATIONS 

We concentrate on those X-ray sources, which demonstrate 
significant thermal emission, own ordinary magnetic fields 
10 n ~ 13 G, own comparatively young ages 10 3-6 yrs, and 
have spins of a few tens of millisecon ds to a few seconds. We 
mainly refer to the collatio n made by Yakovlcv ct al. (2008) 
on cooling neutro n stars, by Becker fc Aschenbachl (|2002l ) on 
X-ray pulsars, bv lHaberll lj2004h on X-ray dim isolated neu - 
tron stars (XDINs), and bv lPavlov. Sanwal fc Teted (|2004l ) 
on central compact objects (CCOs). The sample exhibited 
in Table 1 thus comprises top 17 active pulsar candidates, 
7 XDINs (No. 23-29) and 6 CCOs (No. 5 and 18-22). We 
note that the values of the surface temperatures adopted 
are determined by black body fits, according to the way 
that a SQS emits thermal photons. The XDINs and CCOs 
are considered to be magnetosphere-inactive pulsar candi- 
dates, since their quiescent manifestation on nonthermal ra- 
diation. CCOs, moreover, could only be seen in the soft X- 
ray band without the evidences on the existence of the as- 
sociated PWNe. Additionally, we note that the emission of 
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Figure 2. P — P diagram for the X-ray pulsar sample. The solid 
lines are the death lines both for a typical pulsar (with mass of 
1.4Mq and a radius of 10 km) and for low-mass SQSs (with mass 
of 10 — 4 — 10 — 'Mg). The hexagonal indicate the focused sources 
associated with detectable supernova remnants (SNRs), while the 
pentagons indicate those which exhibit without evident SNRs. 
The triangles, however, indicate 4 XDINs with measurable tempo- 
ral parameters, and the diamonds mar k the CCOs. The timing pa- 
ramet ers for RX J0822-4300 are from lZavlin. Trumper &: Pavlov! 
l|l999h . while the rest are from ATNF Pulsar Catalog 1 . 



the energetic Crab pulsar is much likely to be overwhelmed 
by the nonthermal component originating from its luminous 
plerion, which hampers the detection to the stellar surface 
thermal radiation. As an estimate, we, however, adopt a 
2a upper limit to Crab's surface temperature and an in- 
ferred radius w ith assuming a 2 kpc d istance between Crab 
and the earth (|Weisskopf et alj|2004l ). RX J0822-4300, the 
ce ntral stellar remnant in Puppis A , used to be analyzed 
bv lZavlin. Trumper fc Pavlov! (|l999h basing on the ROSAT 
observations in 1990s. They could not confirm whether the 
X-ray structures surrounding the star belong to the super- 
nova remnant (SNR) or are induced by the probable ac- 
tive magnetosphere. Recent observations on it by Chandra 
and XMM-Newton telescopes, however, did not reveal the 
presence of the associated plerio n, indicating an inert mag- 
netosphere (|Hui fc Becker! 120061 ). Because of the ambiguity 
of this source, we temporarily make it as one of the active 
source candidates. 

Furthermore, we emphasize these X-ray sources that 
own measurable spin parameters in the P — P diagram 
(Fig. 2), by which a distribution of them can be read. 
Other 10 X-ray sources with constraints on the upper lim- 
its of their bolometric luminosity are also denoted in the 
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No. Source t (kyr) T s °5 (MK) (km) T» 2 (MK) ffif (km) L™ ol (10 33 ergs s" 1 ) Refs 



1 PSR B0531+21 1 < 1.85 > 15 ~ 18.8 (1) 
(Crab) 

2 PSR J1811-1925 2 < 1.74 (2) 
(in Gil. 2-0.3) 

3 PSR J0205+6449 0.82-5.4 fa 1.7 fa 2.6 fa 0.44 (3) 
(in 3C 58) 

4 PSR J1119-6127 ~ 1.6 Z-^to.l 3 - 4 ±a3 2 -°tl'.l (4) 
(in G292.2-0.5) 

5 RX J0822-4300 2-5 2.6lt°' 3 g 3.29j;J;^ 5.041°;^ 0.75^ ;^ fa 6.2 (5) 
(in Pup A) 

6 PSR J1357-6429 ~ 7.3 1.7 ±0.2 2.5 ± 0.5 fa 0.37 (6) 

7 RX J0007.0+7303 10-30 < 0.66 fa 12 < 4.0 x 10~ 2 (7) 
(in CTA 1) 

8 PSRB0833-45 11-25 1.06 ± 0.03 5.1±g;| 2.16±g;gf 0.73±g;gf 0-31±g;g| (8) 
(Vela) 

9 PSR B1706-44 ~ 17 2.0lt°'2o 1.81±g;^ fa 0.38 (9) 
(in G343.1-02.3) 

10 PSR B1823-13 ~ 21 1.61±%$ ~ 2.5 fa0.30 (10) 

11 PSR J0538+2817 30 ± 4 2.12±°;°* 1.68 ±0.05 fa 0.46 (11) 
(in S147) 

12 PSR B2334+61 ~ 41 1.62 ± 0.23 1.66l u |g fa 0.11 (12) 
(in G114.3+0.3) 

13 PSR B1916+14 ~ 88 ^to.l °- 8 +5.t ~ °- 03 ( 13 ) 

14 PSRB0656±14 ~ 110 0.91 ± 0.05 fa 14 1.9 ±0.4 fa 0.8 fa 0.96 (14) 
(in Monogem Ring) 

15 PSR J0633±1746 ~ 340 0.5 ± 0.01 8.6 ±1.0 1.9 ±0.3 0.04 ± 0.01 fa 3.2 x 10" 2 (15) 
(Geminga) 

15' 0.482 fa 10 (15') 

16 PSR B1055-52 ~ 540 0.79 ± 0.03 12.31 q% 1-79 ± 0.06 0.46 ± 0.06 fa 0.46 (15) 

17 PSR J2043±2740 ~ 1200 fa 0.9 fa 2 fa 2.0 x 10" 2 (16) 

18 IE 1207.4-5209 ~ 7t_\ 3 1.90 ±0.01 4.5 ±0.1 3.70 ± 0.02 0.83 ± 0.02 fa 2.1 (17,18) 
(in PKS 1209-51/52) 

19 CXOU J232327.9 0.3 6.14 ±0.46 0.41 j;^® i ^ag (19) 
±584842 (in Cas A) 

20 CXOU J085201.4 < 1.1 4.68 ± 0.06 0.28 ± 0.01 0.25 ± 0.02 (20) 
-461753 (in G266.2-1.2) 

21 PSR J1852+0040 5.10 ±3.48 0.9 ± 0.2 3.7 ±0.9 (21) 
(in Kes 79) 

22 PSR J1713-3949 4.4 2.4 15 (22) 
(in G347.3-0.5) 

23 RX J1856.5-3754 ~ 500 fa 0.74 5.0 fa 5.2 x 10" 2 (23) 

24 RX J0720.4-3125 ~ 1300 fa 0.94 fa 6.1 fa 0.21 (24) 

25 RBS 1223 a ~ 400 1.04 0.8 5.1 x 10~ 3 (25) 

26 RX J0420. 0-5022 ~ 110 0.66±°'g4 1.4 2.7xl0 -3 (26) 

27 RX J0806. 4-4123 1.09 ± 6.01 fa0.6 fa 3.6 x 10" 3 (27) 

28 RX J1605.3±3249 1.07 fal.l 1.1 x 10~ 2 (25) 

29 RBS 1774 6 1.04 fal.l 1.1 x 10~ 2 (25) 



1RXS J130848.6±212708 
b 1RXS J214303.7±065419 



Table 1. Observational parameters on X-ray sources, including the surface temperature components T s °^ , 2 , the emission size components 
R^ 2 an d the bolometric X-ray luminosity L^ ol . These are the values detected nearby the earth and fitted or analyzed by the black body 
model. Noting that (i) the age of IE 1207.4-5209 (source No. 18) is adopted according to the estimate of the associated SNR, which is 
given by the reference (17); (ii) the ages of RBS 122 3 and RX J04 2 0. 0-50 22 are assessed by spindown ages in this work; (iii) the age limits 
of sources No. 1, 3-9, 11-12, 14-17, 23-24 are from lYakovlev et al.l l|2008h . while the rest ages are from the corresponding literatures listed 
in the last column; (iv) for 5 XDINs (No. 25 to 29), the black body parameters are obtained with an assumed distance of 100 pc, as a result 
of t he lack or uncertain of thei r distances: (v) th e sp ectral parameters fo r th ese sources are from th e ref erences listed i n the l ast c olumn: 
(1) IWeisskopf et alj i2004 i (2) IZhu et alj <200S|) (3) ISlane et al.l [|2004al) (4) iGonzalez et al.l |2005h (5) iHui fc Becker! d2006l) (6) IZavlin 
1 2007 J > (7l lHalpern et al. d2004^ f8 ~l iMarizali. De Luca fc Caraveol ll2007ll (9) IMcGowan et al.l ll20oi) (lOllPavlov. Kargaltsev fc Brisken 
[120081) (11) IMcGowan et al l <2003) (12) IMcGowan et al.l <l2006h (13) IZhu et alj 20091) (14) IPossenti Mereehetti fc Colpil (119961) (15) 
iDe Luca et al.l ||2005|) (15') Ijackson fc Halpernl ll2005l) (16) IZavlin fc Pavlovl J2004) (17)|Pavlov et alj <2002h (18) |Pe Luca etaT [ 2004 1 ) 
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diagram ((Becker fc Aschenbachl I2002T ) . Besides the death 
line for a typical pulsar, those for low-mass SQSs are 
also drawn in the diagram. These death lines set bound- 
aries of the "graveyard" for pulsars with different mass. 
The timing parameters are from ATNF Pulsar Catalog 
^Manchester et alJ 120051). except for RXJ0822-4300, which 
refers to Zavlin. Triimper fc Pavlovl l|l999l ). 

For magnetosphere-active pulsar candidates, SQSs 
could suggest themselves to reproduce the cooling processes, 
while, for magnetosphere-inactive pulsar candidates, possi- 
ble approaches to understand their current X-ray luminosi- 
ties could also be proposed by SQSs. These two will be pre- 
sented in 63.1 and 53.2. 



3.1 Cooling of active pulsar candidates 

If two isotropic black body emission components — the hot 
component for polar caps and the warm one for the bulk of 
the star — are defined for SQSs, then the equation describing 
cooling processes could be written as (cf. equation ((!])) 



L SH = 47rJ? 2 aT s 4 + 4tt7>T p 4 , 



(11) 



and the relation between T s and T p is given by equation ((6]). 
The luminosity of stellar heating, Lsh, could either follow 
the 1/2-law or the linear-law (see Appendix A). Cooling be- 
havior could thus be calculated by assuming SQSs rotate 
as orthogonal rotators and slow down as a result of mag- 
netic dipole radiation with magnetic field strength at the 
poles of 10 12 G. We provide a parameter space to fit the 
observational temperature-age data, and the comparison be- 
tween the observations and expectations are shown in Figs. 3 
(for 1/2-law case) and 4 (for linear-law case). Some pulsars 
demonstrate two black body components in their thermal 
spectra, e.g. Vela pulsar and the Three Musketeers, implying 
that the temperature inhomogeneity on these pulsars could 
be significant. We thus carried out a temperature-difference 
fit at the same time, as have been shown in the right panels in 
both Figs. 3 and 4. It is worth of noting that, for Geminga, 
if the photon index of the power-law (PL) component is 
thawed when fitting its phase-resolved spectra, the entry of 
the hot black body component could not improve the fits 
remarkably, or it may even becomes an artifact. This may 
mean the fluctuation of the magnetospheric emission during 
a spin cycle might mis lead the understand i ng on its X-ray 
spectra, as analyzed bv ljackson fc Halpernl (|2005l ) and their 
results are denoted by No. 15' in Table 1. If this situation 
holds, the surface temperature fluctuation on Geminga could 
be tiny so that undetectable. We, hence, set a rough upper 
limit to the temperature difference for Geminga of an order 
lower than its surface temperature, as denoted by 15' in the 
"temperature difference- age" figures. 

The proposed relation between the rotational kinetic en- 
ergy loss rate and the bolometric luminosity could inversely 
provide direct measurement to the moments of inertia of 
the active X-ray sources. The results are exhibited in Table 
2 and could be understood in the SQSs' regime, especially 
low-mass ones (cf. the adoption of the coefficients C and rj in 
Figs. 3 and 4). We note that for SQSs with 0.01M Q , O.1M , 
I.OMq and 2.0Mq, the moments of inertia I45 (i.e. moments 



www.atnf.csiro.au/research/pulsar/psrcat. 



Mass 


f (10- 9 M yr- 1 ) 


ibol ( 1Q 38 ergs s -l) 


1O" 4 M 


5.7X1CT 1 


3.5X1CT 2 


1O~ 3 M 


1.2 


7.6xl0" 2 


icr 2 M 


2.7 


1.7X1CT 1 


io- 1 m q 


5.7 


4.7X1CT 1 


1.0M© 


12.4 


2.0 



Table 4. The bolometric X-ray luminosity (Lbol) contributed 
by a SQS's accretion in the propeller regime. The corresponding 
accretion rates, M, are also given. Noting that the values of lumi- 
nosities are calculated by adopting Ae, the latent heat per baryon 
during the phase transition, equals 100 MeV. 



of inertia scaled by 10 45 g cm 2 ) are ~ 2.55 x 10 
~0.55 and ~1.74 respectively. 



-0.01, 



3.2 Thermal X-ray emission of inactive pulsar 
candidates 

As has been discussed in §2.3.2, one probable energy 
source to be responsible for the soft X-ray emission of 
magnetosphere-inactive pulsar candidates lies in the accre- 
tion in the propeller regime. The accretion could be either 
to the interstellar medium or to the fallback disk. Numer- 
ical results for relevant parameters related to this scenario 
will be presented here. Equation {1} for this case could be 
specified as 



GM ■ ?7acc ■ M 

R 



Ae 



T/acc ■ M 



(12) 



with considering equation (j 10|) and the stellar residual ther- 
mal energy is too small to be taken into account. Table 3 
exhibits the parameters for the sources under the scenario 
of accretion in the propeller regime. These parameters in- 
clude the light cylinder radius m, the corotation radius r co , 
and the Alfven radius r m . Table 4 lists the luminosities of 
SQSs under the accretion scenario. It could be concluded 
that the sources' current observational thermal X-ray lumi- 
nosities could be interpreted by considering the parameters 
H, the accretion rate in unit of that of Eddington, and ?7acc, 
the accretion efficiency, e.g. fj, = 0.01, r^cc = 0.001. Noting 
that for the enigmatic central stellar remnant, RX J0822- 
4300, we also present its results in this case, since we could 
not exclude the possibility that the star's bolometric lumi- 
nosity is of accretion origin. 



4 CONCLUSIONS AND DISCUSSIONS 

We collate the thermal observations of 29 X-ray isolated pul- 
sars and, in the SQSs' regime, for the magnetospherically 
active pulsar candidates, establish their cooling processes 
(Figs. 3 and 4), while for the magnetospherically inactive or 
dead pulsar candidates, interpret the X-ray luminosities un- 
der the accretion scenario (Tables 3 and 4) . SQSs, because of 
the possibility of being low-mass, could provide an approach 
to understand the small black body emission sizes. We note 
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Figure 3. Left panel: Cooling curves for SQSs, if the 1/2-law between the bolometric luminosity and the spin energy loss rate holds. Right 
panel: Corresponding temperature differences between the hot and warm components of SQSs, or AT = — T s °°. The parameters 
in both panels: M = O.IMq, C = 10 16 (solid lines); M = I.OMq, C = 10 16 (dashed lines); M = I.OMq, C = 10 15 (dotted lines); 
M = O.OIM0, C = 10 15 (dash-dot lines) (C is in unit of ergs 1 / 2 s -1 / 2 ). For two curves with same M and C, the upper one corresponds 
to an initial spin of 10 ms, while the lower one 100 ms. Noting that the errors on the surface temperatures of PSRs J0205+6449 (No. 3) 
and J2043+2740 (No. 17) are not provided by the authors of the references, we then conservatively adopt them as deviating from the 
central values by a factor of 2. 
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Figure 4. Left panel: Cooling curves for SQSs, if the linear-law holds. Right panel: Temperature differences for this case. The parameters: 
M = I.OMq, T) = 0.01 (solid lines); M = 1.0M Q , 77 = 0.001 (dashed lines); M = O.IMq, 77 = 0.1 (dash-dot lines); M = 0.01M©, r? = 0.1 
(dotted lines). As in Fig. 3, for two curves with same M and rj, the upper one corresponds to an initial spin of 10 ms, while the lower 
one 100 ms. 



that for SQSs with mass of 0.01M Q , 0.1M o and 1M , their 
radii are ~1.8, ~3.8 and ~8.3 km respectively. On the other 
hand, a linkage between the pulsars' rates of rotational ki- 
netic energy loss and the bolometric X-ray luminosities is 
explored by SQSs (see Appendix A), and Table 2 exhibits 
the evaluated moments of inertia. We hence conclude that 
SQSs could not be ruled out by the observations on the X- 
ray isolated pulsars. Additionally, the various performances 
of X-ray pulsars may indicate their current states or proper- 
ties and imply their possible evolutionary history. We thus 
extend a discussion below. 



4.1 Spin-powered pulsars 

Spin has always been a significant energy source for ac- 
tive pulsars, by which multiwave bands' nonthermal radi- 



ation are driven, including those both of the pulsars' and 
the surrounding plerions'. Some phenomenological studies 
demonstrate certain regularities for such a population. A 
brief summary is that the nonthermal X-ray luminosity 
and spin energy loss rate own a relation of Lx = W~ 3 E 
(|Becker fc Triimper|[l997l ). and the 7-ray lumin osity is pro- 
porti onal to the square root of E or L 7 oc E 1/2 l|Thompsonl 
Il997l ). We, additionally, note that V oc E 1/2 , where V is 
the potential drop along the open field region. Therefore, it 
seems that for such a population the more younger the larger 
E and V that a pulsar owns (cf. Fig. Al, bottom panel) and 
the more luminous of the radiation in the hard X-ray and 
7-ray bands. 

Besides nonthermal emission, spin may also be an en- 
ergy origin of pulsars' thermal radiation as the relations of 
L™ ol oc E 1/2 or L™ ol ~ lO" 3 ^ could exist observationally 
(see Appendix A for details). SQSs predominantly follow 
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No. Source v (s _1 ) v (s -2 ) 7J 5 ■ C 2 (ergs s _1 ) J 4e 



1 


PSR B0531+21 (Crab) 


30.225437 


-3.862X10 -10 


7.7xl0 29 


4.1x10" 


-5 


2 


PSR J1811-1925 


15.463838 


-1.052xl0 -11 








3 


PSR J0205+6449 


15.223856 


_4.495xl0 -11 


5.8xl0 28 


4.6x10" 


-5 


4 


PSR J1119-6127 


2.452508 


-2.419X10" 11 


1.7xl0 30 


8.5x10" 


-4 


5 


RX J0822-4300 


13.2856716499(3) 


-2.6317(3) xlO" 11 


2.7xl0 30 


2.8x10" 


-4 


6 


PSR J1357-6429 


6.020168 


-1. 305x10-" 


4.4xl0 28 


1.2x10" 


-4 


7 


RX J0007.0+7303 


3.165922 


-3.623X10" 12 


3.5xl0 27 


8.8x10" 


-5 


8 


PSR B0833-45 (Vela) 


11.194650 


-1. 567X10" 11 


1.4xl0 28 


4.5x10" 


-5 


9 


PSR B1706-44 


9.759978 


-8.857xl0~ 12 


4.2xl0 28 


1.1x10" 


-4 


10 


PSR B1823-13 


9.855532 


-7.291X10" 12 


3.2xl0 28 


1.1x10" 


-4 


11 


PSR J0538+2817 


6.985276 


-1.790X10" 13 


4.3xl0 30 


9.3x10" 


-3 


12 


PSR B2334+61 


2.018977 


-7.816XKT 13 


1.8xl0 29 


1.7x10" 


-3 


13 


PSR B1916+14 


0.846723 


-1.523XHT 13 


1.8xl0 29 


5.9x10" 


-3 


14 


PSR B0656+14 


2.598137 


-3.713X10" 13 


l.lxlO 31 


1.7x10" 


-2 


15 


PSR J0633+1746 (Geminga) 


4.217640 


-1.952X10" 13 


3.2xl0 28 


9.9x10" 


-4 


16 


PSR B1055-52 


5.073371 


-1.501X10- 13 


6.9xl0 30 


1.5x10" 


-2 


17 


PSR J2043+2740 


10.402519 


-1.374xl0- 13 


7.1X10 27 


3.5x10" 


-4 



a J 45 = //(10 45 g cm 2 ) 

Table 2. The moments of inertia of active pulsar candidates suggested by SQSs. The 5th column lists the values derived by the 1/2-law, 
while the 6th column lists those obtained by the linear-law. The temporal parameters are fr om ATNF Pulsar Catalog ( s ee foo tnote 1 for 
the website). For RX J0822-4 300, we note that we adopt the spin parameters obtained bv lZavlin. Triimper S z Pavlov (1999) basing on 
a 4.5 yr-span data set, though [Hui fc Becker (20o3) provided more recent values using a 0.5 yr-span data set. 



these relations to accomplish their cooling processes, since 
their residual thermal energy could be quite inadequate to 
sustain a long term X-ray thermal emission. This could be 
the different nature when comparing with neutron stars, 
whose cooling is powered by residual thermal energy. 

Furthermore, it could be an intrinsic distinction for 
SQSs from neutron stars that the heat capacity of SQSs is 
much lower than that of neutron stars. The thermal X-ray 
emission of SQSs could only be sustain ed by heating, a nd the 
missing compact object in SN1987A (|McCrav| 12007) could 
be a SQS with negligible residual thermal energy and thus 
is not X-ray loud observationally. 



4.2 Accretion-powered pulsars 

XDINs and CCOs could be representative populations of 
magnetosphere-inactive pulsars. Nevertheless, their energy 
origin that could power the X-ray emission is still an enigma. 
The bolometric luminosities of IE 1207.4-5209, RX J1856.5- 
3754 and RX J0720. 4-3125 exceed their spin energy loss 
rates by a factor of ~ 60, ~ 15 and ~ 45, respectively; 
they do not appear to be powered by spin. 

CCOs could be a group o f weakly-mag n etized 
and long-initial-period pulsars (|Pavlov et al.l |2002| ; 
iGotthelf fc Halpern|[2007l ), because of which the potential 
drop along the open field lines could be much less than 
~ 10 12 V, so that the primary plasma could not be fully 
accelerated and result in an inactive magnetosphere. The 
black body fits to the thermal spect ra of CCOs always result 
in ex tremely small emission sizes (|Pavlov. Sanwal fc Teterl 
|2004| ). implying possible existence of low-mass SQSs, even 
as low as ~ KT 4 M Q (|Xul l2005h . Hence, from another 
perspective, CCOs' E might be very low, so that effective 



magnetospheric emission and plerions could not be driven. 
CCOs could thus be the representatives of a group of natal 
inactive pulsars. Accretion to the fallback ejecta of the 
associated SNRs could power their soft X-ray luminosity, 
and th is scenario have be e n men t ioned by several authors 
(e.g., Chakrab artv et ail l200l|; [k argalts ev et al. 2002; 



iFesen. Pavlov fc SanwallbOOfj ; lHalpern et al.ll2007l ) 

XDINs, owning comparatively larger spindown ages and 
stronger magnetic fields, might be the descendants of mag- 
netars. Considering sim ilar properties between them l|Haberll 
12004 lMereghettill2008l ). one could not exclude that XDINs 
are still being powered by the decaying magnetic fields. How- 
ever, observations have revealed the probable existence of 
residual disks around such p ulsars, and thus their radiation 
may be of accretion origin (|Lo Curto et alj|2007l ). In this 
case, XDINs could be the evolved products of active pul- 
sars, and thanks to the accretion so that they are still visi- 
ble after their death. Most XDINs have absorption features 
in their spectra, being similar to those of CCOs. Hence, if 
they are really undergoing accretion in the propeller regime, 
their spindown rates, P, could then be of accretion ori- 
gin rather than causing by the magnetic dipole braking. 
Thus their magnetic fields could be much lower than the 
values determined by canonical magnetic dipole radiation, 
since the absorption features could be electron cyclotron 
lines. Therefore, the common properties of absorption fea- 
tures between XDINs and CCOs could imply that XDINs 
might be older CCOs. XDINs and CCOs have not mani- 
fested themselves as radio pulsars, their radio-quiet demon- 
strations could be an intrinsic property rather than being of 
beaming origin, since they might be dead or natal inactive 
ones. So they could be grouped as radio-quiescent pulsars. 
For RX J0007. 0+7303, an active pulsar candidate, its lack of 
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No. Source v (s^ 1 ) v (s -2 ) r L (10 10 cm) r co (10 8 cm) r m ■ /i 2 / 7 (10 8 cm) 



5 


RX J0822-4300 


13.2856716499(3) 


-2.6317(3)xl0~ n 


3.59xl0~ 2 


1.24xl0~ 2 


5.05xl0 _1 












2.67xl0~ 2 


8.74xl0 _1 












5.76xl0~ 2 


1.51 












1.24xl0 _1 


2.62 












2.67xl0~ 1 


4.53 


18 


IE 1207.4-5209 


2.357764 


(-3.669 ± 5.1) x 10~ 16 


0.20 


0.18 


0.27 












0.39 


0.47 












0.85 


0.82 


19 


CXOU J232327.9+584842 












20 


CXOU J085201. 4-461753 












21 


PSR J1852+0040 


9.531742 


(3.089 ± 2.5) x 10~ 14 








22 


PSR J1713-3949 


2.548086 










23 


RX 11856 5-3754 


141739 


—5 980 x 10 -16 


3.37 


1.19 


3.49 












2.56 


6.04 












5.52 


10.46 


24 


RX J0720.4-3125 


0.119174 


-9.918 x 10" 16 


4.00 


1.33 


4.68 












2.88 


8.10 












6.20 


14.02 


25 


RBS 1223 


0.096969 


-1.053 x 10~ 15 


4.92 


1.53 


4.97 












3.30 


8.59 












7.11 


14.86 


26 


RX J0420.0-5022 


0.289597 


-4.193 x 10~ 14 


1.65 


0.74 


6.38 








±8.0 x 10~ 13 




1.59 


11.03 












3.43 


19.08 


27 


RX J0806. 4-4123 


0.087943 


(3.558 ± 5.0) x 10~ 14 








28 


RX J1605. 3+3249 


0.1453? 










29 


RBS 1774 


0.105966 











Table 3. Numerical values of the parameters under the propeller regime if the sources own detectable spindown parameters. The 
parameters include the light cylinder radius r^, the corotation radius r co , and the Alfven or the magnetosphere radius r m . For each 
source, three values on r co and r m are given, which, from top to bottom, correspond to the results of a O.OIMq/O.IMq/I.OMq-SQS 
respectively, but the results of RX J0822-4300 are presented more d e tailed star ting from the case of a 1Q- 4 M W - SQS. The temporal 
parameters for RX J1605.4-4123 and RX J0822-4300 refer to lHaberll d2004) and IZavlin. Trumper fc Pavlovl l ll999h respectively, while 
those of the rest are from ATNF Pulsar Catalog (see footnote 1 for the website), and we note the errors when they are significant enough. 



radio signal could be the result of an unfavorable geometry, 
or its radio beam sweeps away f r om th e earth as has been 
analyzed bv lBrazier fc Johnston! |l999ri . 



The similar thermal manifestations of accreting pul- 
sars (e.g. XDINs and CCOs) and cooling pulsars may have 
caused a confusion about distinguishing between the two 
classes (|Treves et al.ll2000l) . However, the activeness of the 
magnetosphere could provide a way to achieve an iden- 
tification. If pulsars are in fact SQSs with rapid rota- 
tion, cooling pulsars are then likely to be undergoing the 
spindown-powered heating evolution, and multi-bands non- 
thermal emission originating from luminous magnetospheres 
and even PWNe will accompany such coolers during the pro- 
cesses. In contrast, accreting pulsars may have not demon- 
strate such characteristics observationally. In this context, 
the CCO in the SNR Cassiopeia A (source No. 19 in Table 
1) should be an accretor because of the lack of PWN a nd 
nonthermal power-law component (jPavlov fc Lu na 2009). 



APPENDIX A: SPIN-POWERED THERMAL 
EMISSION? 

In the solid quark star regime, magnetosphere-active pul- 
sars' cooling depends predominantly on the stellar heating 
as a result of the lack of sufficient residual heat as has been 
analyzed in §2. The pulsar activity induced stellar heating 
would have intrinsically set up a linkage between the X-ray 
bolometric luminosity and the spin energy loss. We, hence, 
present here a phenomenological study on such a relation. 

The sample of active p ulsars firstly includes the top 
17 sources listed in Table 1. iBecker fc Aschenbachl (2002) 
summarized X-ray pulsars, in which other 10 sources' upper 
limits on the bolometric luminosities could be defined. These 
10 pulsars are meanwhile considered. Fig. Al (top panel) 
illustrates these sources' bolometric luminosities L^ ol as a 
function of their spin energy loss rates E, and Table Al lists 
the fitted parameters led by the function of 

logi^ ol =pl-log£ + p2. (Al) 

The fits are carried out both only for the top 17 pulsars taken 
from Table 1 (as Group A) and for the whole sample with 
the other 10 sources included at the same time(as Group 
B). The best fits, moderate clearly, suggest a 1/2- law on the 
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Group pl a p2 Corr. Coef. 6 \l (d.o.f) c 



A. 0.4561 ±0.2315 16.20 ± 8.30 0.3277(14) 

(1) -3.280to494 0.7487 0.8607(15) 

B. 0.5918t°jg« n.66±™J°° 0.6081(24) 

(1) -2.896 ±0.403 0.7730 0.9964(25) 



a Values in the parenthesis are frozen during the fits. The errors are in 95% confidence level. 
b Correlation coefficient of the data between bolometric luminosity and E. 
c Reduced x 2 > or \ 2 P er degree of freedom (d.o.f). 

Table Al. Fitting parameters for L^ ol — E data, which has been shown in Fig. Al (top panel). 



relation between pulsars' bolometric luminosities and the 
spin energy loss rates, or 

LZi{E) = CE 1/2 , (A2) 

with a coefficient C = 10 p2 in unit of ergs 1 / 2 s -1 / 2 . If pi is 
fixed at 1, a linear-law then is implied, i.e. 



LZi{E) =r,-E, 



(A3) 



with a coefficient rj = 10 p2 , or, in other words, r\ is the 
conversion efficiency. The fits obtain r\ 



10" 



which is 



simil ar to the nonthermal X-ray case (|Becker fc Trumperl 
1 19971 ). For the 1/2-law case, the conversion efficiency turns 
out to be the function of E, and it could be 



V(E) 



CE 1 ' 2 
E 



C 

£1/2' 



(A4) 



which could in turn, besides the fits, provide a reference 
to the value of the coefficient C by a natural constraint of 
r](E) < 1. Taking B0823±26, the minimum E in our sam- 
ple, as an example, its C could be less than ~ \/l0 32 = 10 16 
ergs 1 / 2 s~ 1,/2 . Additionally, we note here that we also present 
a age-i? relation in the bottom panel of Fig. Al, which illus- 
trates a natural trend for spin-powered pulsars. 
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Figure Al. Functions L^AE) (top panel) and t(E) (bottom 
panel) of active pulsar candidates. In the top panel, the fits are 
carried out both for Group A (including top 17 pulsars listed 
in Table 1; they are marked by dark points and their numbers 
here) and Group B (including the members in Group A and the 
other 10 sources whose upper limits on the bolometric luminosi- 
ties could be defined observationally; the 10 sources are marked 
by g re y points and letters). Red lines are the fitted results for 
Group A, while blue lines are those for Group B. In both groups, 
solid lines provide the best fit to the data, while the dashed lines 
give the fits by freezing pi at 1. We note that the 10 sources with 
uppe r limits on their bolometric luminosities defined are taken 
from [Becker fc Aschenbachl {2002), and they are a. B1509-58, b. 
B1951+32, c. B1046-58, d. B1259-63, e. B1800-21, f. B1929+10, 
g. B0540-69, h. B0950+08, i. B0355+54, j. B0823+26. 
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